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First-principles calculations are used to calculate the strain dependencies of the binding and 
diffusion-activation energies for Ge adatoms on both Si(001) and Ge(001) surfaces. Our calculations 
reveal that the binding and activation energies on a strained Ge(001) surface increase and decrease, 
respectively, by 0.21 eV and 0.12 eV per percent compressive strain. For a growth temperature 
of 600 °C, these strain-dependencies give rise to a 16-fold increase in adatom density and a 5-fold 
decrease in adatom diffusivity in the region of compressive strain surrounding a Ge island with a 
characteristic size of 10 nm. 



In heteroepitaxial growth of lattice-mismatched thin 
films, the Stranski-Krastanov (SK) growth mode has 
been widely investigated as a basis for self-assembling 
arrays of coherent nanostructured islands, commonly re- 
ferred to as quantum dots (QDs)ii The need for highly 
monodisperse QD arrays in semiconductor optoelectronic 
device applications has motivated extensive research into 
the microscopic mechanisms influencing the evolution of 
island size distributions during SK growth. While the 
stabilizing thermodynamic effects associated with the 
elastic interactions between strained islands have been 
investigated in some detail (e.g., Refs. HI and refer- 
ences cited therein), the role of various proposed kinetic 
mechanisms on the development of island size distribu- 
tions remains less clear. 

Within a self-consistent mean-field rate theory, Ko- 
duvely and ZangwillA demonstrated that with decreas- 
ing island-island separation, a strain-mediated decrease 
in the barrier for adatom-island detachment leads to a 
reduction in the mean island size, with an associated 
narrowing of the size distribution. These findings are 
qualitatively consistent with experimental observations 5 
in InAs/GaAs. Madhukar 6 and Penev, Kratzer and 
SchcfheBi^ have considered the growth of InAs QDs on 
GaAs, where increasing island size leads to a buildup of 
compressive elastic strains in the surrounding substrate £ 
Within simplified models of diffusion-limited growth 
these authors demonstrated that, in the InAs/GaAs het- 
eroepitaxial system, the strain-dependence of the param- 
eters governing adatom diffusion gives rise to a reduction 
in the flux of adatoms reaching larger islands relative to 
small ones, leading to a reduced rate of coarsening and 
an associated narrowing of the size distribution. 

While the potentially important consequences for is- 
land growth kinetics arising from strain dependencies 
in adatom binding and migration energies have been 
clearly demonstrated, attempts to determine the mag- 
nitude of these effects in specific systems have been un- 
dertaken in relatively few systemsi&iSiiiiiSii&ii These ef- 
fects are most readily investigated using first-principles 
calculationa&iSiiAii 5 ^^, as they are extremely difficult to 
isolate experimentally. The purpose of the present work 
is to investigate the effect of strain-dependent adatom 
diffusion and binding energies upon island growth kinet- 



ics in Ge/Si(001). This system represents one of the most 
widely studied examples of QD formation induced by 
SK growth; 17 i 18 i 19 i 20 i 2:L i 22 yet to date the effect of strain 
upon Ge adatom binding energies and diffusion rates on 
Ge wetting layers remains unstudied. We employ first- 
principles calculations to compute the strain dependence 
of Ge adatom energetics on Si(001) and Ge(001) surfaces. 

Our calculations were set up as follows. We con- 
sider adatom diffusion in the dilute limit (as opposed 
to dimer diffusion 2 ^). The initial atomic configuration 
for all of our calculations is taken to be the well-known 
minimum energy c(4 x 2) reconstruction (depicted in the 
upper left corner of Fig. The adatom calculations 
are performed using a supercell geometry where an ar- 
tificial three-dimensional periodicity is imposed on the 
system. This artificial periodicity is harmless, provided 
that convergence with respect to the distance between 
the periodic images is achieved. Following earlier studies 
of adatom diffusion of similar systems j^S 5 *^ our super- 
cell consists of two repetitions of the surface unit cell, 
depicted in Fig. ^ along the dimer rows direction (i.e. 
a 4 x 4 supercell). Twelve layers of atoms were used, 
separated by 6 layers of vacuum. By varying the thick- 
ness of the slab and the width of the vacuum separating 
the periodic images perpendicular to the surface, it was 
verified that the selected system size provides energies 
with an accuracy of the order of 0.03 eV. Single adatoms 
were placed on each of the two slab surfaces, which has 
the effect of doubling the adatom's energetic contribu- 
tion, thereby improving the accuracy of the calculations. 
All calculations were performed using the ab initio total- 
energy and molecular-dynamics program VASP (Vienna 
ab initio simulation package) developed at the Institut 
fur Materialphysik of the Universitat Wien2L2&, which 
implements Vanderbilt ultra-soft 29 pseudo-potentials 30 . 
The energy cutoff for the plane-wave basis set was set to 
150 eV and a 2 x 2 x 1 mesh of fc-points was used for the 
4x4 supercell. All atoms were allowed to relax. 

On the Si (001) surface, the diffusion of a Ge adatom 
along the dimer rows is known to be of the order of 1000 
times faster than across dimer rows^i and we therefore 
focus solely on the diffusion along dimer rows. Since a 
typical surface consists of terraces where the direction of 
the dimer rows changes by 90° at each monoatomic step, 
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FIG. 1: Geometries of the binding site and saddle point con- 
figurations for a Ge adatom (filled circle) diffusing on the Si 
(001) and Ge (001) surfaces. Only the two topmost atomic 
layers are shown. The size of each circle reflects the atom's 
proximity to the observer. In the lower left quadrant, crosses 
indicate the adatom position for the 6 other candidate binding 
sites that were considered. 



fast diffusion in any direction is possible at the meso- 
scopic level, even though the fast diffusion is unidirec- 
tional at the microscopic level. Earlier studies^2£i2£i22 
have unambiguously determined the location of the bind- 
ing site and the activated state of the Ge adatom on the 
Si (001) surface (see Fig. as well as the precise configu- 
ration of the surface dimers in the vicinity of the adatom. 

Since experimental and computational evidence is 
scarcer in the case of the Ge adatom on the Ge (001) 
surface,? 3 we considered 7 possible binding sites and 
found the minimum energy site to be as shown in Fig- 
ures Q] and This site remains energetically favored for 
values of the substrate lattice parameter up to 3% larger 
than the Si lattice parameter. Beyond that threshold, 
a binding site analogous to the binding site of Ge on Si 
(001) becomes favorable. We shall neglect this alternate 
binding site, as the substrate lattice parameter required 
to stabilize it falls outside the range sampled during Ge 
on Si (001) heteroepitaxy. 

We identified the diffusion path using the nudged elas- 
tic band metho d 34 ^ and refined the position of the sad- 
dle point (see Figures ^ an( l EJ using the quasi-Newton 
algorithm^ We focused on the diffusion along the val- 
leys between the dimer rows, because the binding sites 
located atop the dimer row had a binding energy at least 
0.3 eV larger than the saddle point energy for diffusion 
within the valleys, strongly suggesting that any other 
saddle points located on the dimer rows would also have 
a higher energy than the saddle point we analyzed. The 
nearly onc-dimensional diffusion and the adatom's pref- 
erence for sites located in the valley between the dimer 
rows agrees qualitatively with the results of earlier cal- 
culations based on semi-empirical potentials^ although 
the precise location of the binding sites and of the saddle 
points differ. 

For each type of surface (Si or Ge), the energies of 
three geometries were calculated (the free surface, the 
saddle point and the binding site configurations) at var- 
ious levels of biaxial strain imposed parallel to the plane 




FIG. 2: Geometries of the binding site (left) and saddle point 
(right) configurations for a Ge adatom (indicated by an arrow) 
diffusing on the Ge (001) surfaces. 
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FIG. 3: Strain-dependence of the binding energies Et (left), 
the saddle point energies E a + Eb (middle) and the activation 
barriers E a (right) for a Ge adatom diffusing on Si(001) (top) 
and Ge(001) (bottom) surfaces. These plots are obtained by a 
polynomial fit to the calculated energies as a function strain. 
Let a, asi and og c respectively denote the substrate, the Si, 
and the Ge lattice parameters. 



of the surface. The resulting energy versus strain re- 
lationships are plotted in Fig. [3] Our results reveal 
four important observations. First, linear approxima- 
tions to the strain-dependence of the binding energies 
and activation barriersii must clearly be used with care. 
Second, the Ge adatom binding energy does not nec- 
essarily exhibit a minimum when the lattice parameter 
of the substrate matches the lattice parameter of bulk 
Ge. Third, the sign of the change of these energies 
under strain are highly system-dependent: calculations 
on the structurally similar In/GaAs(001) heteroepitaxial 
systemS. found the strain-dependences of the binding (Eb) 
and saddle-point energies (Eb + E a ) the be of a sign op- 
posite to the ones in the Ge/Ge(001) system. The sign of 
the strain-dependence of the activation barrier (E a ) for 
Ge on both Si (001) and Ge (001), however, agrees with 
earlier studies in In/GaAs(001)£ and Si/Si(001) 11 T 12 i 13 i 23 . 
Finally, the magnitude of the effect of strain on bind- 
ing and saddle point energies found in the case of Ge on 
Ge(001) is by far the largest, relative to any other system 
studied so far, 8 . 10 i n . 12 i 13 i 23 

In order to quantify the importance of these results in 
the context of quantum dot growth, we calculated the 
strain field in the vicinity of a Ge (105)-terminated pyra- 
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FIG. 4: Strain field experienced by the substrate surface in 
the vicinity of a Ge island (dotted lines) of a typical size (10 
nm) terminated by (105) facets. The homogenous component 
of the in-plane strain relative to the Si lattice parameter is 
plotted ((e xx + e yy )/2, for a surface normal to the z axis). 



midal island with a characteristic size of 10 nmS For 
the calculation of strain fields around coherent strained 
islands, both finite-element techniques (e.g., Ref. 0) as 
well as approximate analytical methods^ have been ap- 
plied previously. In the present work, elastic strain fields 
have been derived from the results of the linear stabil- 
ity analysis of Spencer, Voorhees and Davisi^ In this 
analysis, the linearized strain fields arising from Fourier- 
mode perturbations in the surface height were derived 
within isotropic elasticity theory. By summing the re- 
sulting expressions over the Fourier amplitudes describ- 
ing the shape function of a Ge pyramid, we calculated 
the in-plane strain field (shown in Fig. at the surface 
of a Ge wetting layer (three atomic-planes in thickness) 
on Si(001). In Fig. 0] the strain is referenced to the Si 
lattice constant, so that zero corresponds to a wetting 
layer epitaxially strained on the Si substrate. In agree- 
ment with previous calculations for related systemSf&Si^l 
the strain field surrounding the island is compressive in 
nature. 

In Fig. the strain is seen to range in magni- 
tude from zero to roughly —1% in the vicinity of the 
Ge(105)-terminated island. Over this range of compres- 
sive strains, the binding energy of a Ge adatom on a 
Ge surface varies by about 0.21 eV, while the activa- 
tion barrier varies by about 0.12 eV, as seen in Fig. [3] 
At a typical deposition temperature of 600 °C, these 
variations correspond, respectively, to a 16-fold varia- 
tion in equilibrium adatom density and a 5-fold varia- 
tion in adatom mobility. In contrast, the correspond- 
ing strain-dependences on the Si (001) surface are much 
less pronounced. Of course, the strain-dependence of 
the entropic prefactors^ could very well modify those 
order-of-magnitude estimates based solely on the strain- 
dependence of the energetic contributions. 

To further explore the implications of our calculated 
results, we follow the analysis of Penev, Kratzer and 



SchefilerySi employing a simple surface diffusion-limited 
one-dimensional model of adatom diffusion between two 
neighboring islands. We expand slightly on their work by 
including boundary conditions at the island edges that 
impose local equilibrium between adatoms and islands. 
For simplicity, we consider all entropic prefactors to be 
strain- independent. We perform a stability analysis by 
considering two islands of equal size located at x = 
and x = I and calculate the changes in the adatom flux 
toward each island as the size of one is perturbed. Let 
F ,Fi denote the flux towards each island and consider 
the (dimensionless) asymmetry in the flux towards each 
island, F = — Fo) / (2(j>l), where <j) is the adatom de- 
position rate and I is the distance separating the two 
islands. Solving the above diffusion problem leads to 

K (e^M-e^Q) Mi 

M (j>l I M ' U 

where 

M n = -L J* ^_ L^ n eP (E a{x)+ E b[x))dx _ (2) 

and where E a (x) and Ef, (x) are, respectively, the activa- 
tion barrier and binding energy as a function of position 
x; Ei.o and E^i are, respectively, the energies of an atom 
bound to the islands located at x = and x = I; K is 
a constant incorporating all entropic prefactors and /3 is 
reciprocal temperature (ksT) . 

Our analysis focuses on the diffusion of Ge adatoms 
on the Ge (001) surface, strained to match the lattice 
constant of the Si substrate, since a Ge wetting layer 
covering the Si substrate is known to form in SK growth. 
We assume that adatoms are insensitive to the presence 
of Si under the Ge wetting layer, that no substantial Si- 
Ge interdiffusion occurs and that surface reconstructions 
are unaffected by epitaxial strain. Accounting for the 
latter (e.g., through the inclusion of a strain-dependent 
concentration of "missing dimers"— ) may provide an- 
other source of strain-dependent diffusion deserving fur- 
ther consideration. 

Under the above assumptions, the quantity E a (x) + 
Eb (x) decreases under a compressive strain, as shown 
in the middle panels of Fig. [3] When the two islands 
are identical, F = 0, since Eij — Ei$ and M\ = (as 
E a (x)+Eb (x) is then symmetric with respect to x = |). 
Now, consider an increase in the size of the island at 
x = I. The first term of Equation JTJ describes the ther- 
modynamic driving force for coarsening: e° Ei < 1 decreases 
due to capillarity (a larger island has a smaller surface to 
volume ratio). This term thus causes F to become pos- 
itive, increasing the flux towards the larger island. This 
term is inversely proportional to cf>, demonstrating that 
increasing the deposition rate would act to reduce this 
natural coarsening effect. 

The second term of Equation which is indepen- 
dent of deposition flux, quantifies the effect of strain- 
dependent adatom binding and migration energies. A 
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slight increase in the size of the island at x = I induces a 
compressive strain in the wetting layer in the vicinity of 
that island, thus inducing a decrease in the saddle point 
energy E a (x) + Eb(x) in the same area and causing Mi to 
become negative. The result is a relative increase in the 
adatom flux towards the larger island. In other words, 
this analysis suggests that the strain-dependencies of Ge 
adatom binding and migration energies plotted in Fig. [3] 
act to accelerate the rate of coarsening of larger islands. 

In summary, first-principles calculations have been em- 
ployed to compute the strain dependence of Ge adatom 
binding energies and activation energies for diffusion on 
both Si(OOl) and Ge(OOl) surfaces. For a growth tem- 
perature of 600 °C, these strain-dependencies give rise 
to a 16-fold increase in adatom density and a 5-fold de- 
crease in adatom diffusivity in the region of compressive 
strain surrounding a Ge island with a characteristic size 
of 10 nm. Within a simplified model of diffusion-limited 
growth, these strain dependencies are found to have the 



qualitative effect of accelerating the natural coarsening 
rate of larger islands, in contrast with earlier results ob- 
tained for the related InAs on GaAs system^ The large 
magnitude of the strain dependence of adatom energetics 
on Gc(001) obtained in the present calculations also con- 
trasts with earlier findings in other systems&iS*ii*i2i±2i2£ 
and warrants further consideration of these effects in 
more detailed kinetic models to further elucidate their 
effect on island growth. 
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